Parametric level statistics in random matrix theory: Exact solution 
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An exact solution to the problem of parametric level statistics in non-Gaussian ensembles of 
N X N Hermitian random matrices with either soft or strong level confinement is formulated within 
the framework of the orthogonal polynomial technique. Being applied to random matrices with 
strong level confinement, the solution obtained leads to emergence of a new connection relation 
that makes a link between the parametric level statistics and the scalar two-point kernel in the 
thermodynamic limit. 
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In recent years, the theory of non-Gaussian random 
matrix ensembles has experienced a sound progress mo- 
tivated by new ideas in quantum chromodynamics 
(QCD) and mesoscopic physics. Both fields, concentrat- 
ing on physically different objects which besides have 
incommensurable spatial scales, make use of the same 
language of random matrix theory j^] to reveal the uni- 
versal statistical manifestations of symmetries inherent 
in a physical system. In QCD, invariant non-Gaussian 
random matrix models incorporating global symmetries 
of the QCD Dirac operator appear to describe the uni- 
versal features of low-energy part of the spectrum 
of the Dirac equation. In mesoscopic physics, similar 
matrix models of appropriate symmetries arise in iden- 
tifying the universal electron level statistics in normal- 
superconducting hybrid structures [Q, as well as in the 
context of mesoscopic electron transport Q . A new class 
of non-Gaussian random matrices |^ with finite level 
compressibility should also be mentioned; these serve as a 
basis for constructing the toy models of universal spec- 
tral statistics expected at the edge of the metal-insulator 
transition. 

One of the most promising ways of studying the non- 
Gaussian random matrix ensembles consists of revealing 
the universal connection relations (i) between different 
entities (statistics) within the same matrix model, and 
(ii) between the same statistics in different matrix mod- 
els (which only differ by underlying symmetries). In both 
cases the connection identities (if they exist) originate 
from underlying mathematical structure of random ma- 
trices which deserves to be studied in its own right, also. 
Being brought together, these two approaches have led to 
a substantial breakthrough in understanding the univer- 
sality phenomenon in invariant random matrix ensembles 
characterized by strong level confinement Q . 

(i) The Shohat method of Refs. §[10 1 is a typical ex- 
ample of the first approach. This formalism, employ- 
ing the free fermion representation of the joint distri- 
bution function of eigenvalues in unitary invariant ran- 
dom matrix ensembles, reduces the problem of finding 
the n-point spectral correlators in the thermodynamic 
limit to a simpler problem of solving an effective one- 
particle Schrodinger equation. In the absence of a fine 



tuning of confinement potential to so-called multicritical 
points pT| , this equation involves two ingredients: the 
macroscopic mean density of states vi^{E) (which is a 
one-point spectral characteristics) and the term deter- 
mined by the singularities of confinement potential. In 
this way, one relates the n-point spectral characteristics 
to the one-point characteristics. 

(ii) Recent studies by §ener and Verbaarschot [ p^ , 
and by Widom ||l^ are representatives of the second 
approach. These authors, by using the two completely 
different constructions (of skew |lj] and ordinary or- 
thogonal polynomials, respectively) succeeded in deriving 
universal relations between eigenvalue correlations in a 
class of non-Gaussian random matrix ensembles with the 
same, finite polynomial, confinement potential but with 
different symmetries. These remarkable relations express 
the spectral correlation functions of matrix models with 
orthogonal {(3 — 1) and symplectic (/3 = 4) symmetries in 
terms of the scalar two-point kernel of the corresponding 
ensemble with unitary symmetry [(3 = 2). 

Contrary to conventional spectral statistics, which has 
received a detailed study in the literature, quite a bit is 
known | [l^ , [r^ about parametric level statistics in non- 
Gaussian random matrix models. In the present Letter, 
we address this problem for ensembles with unitary sym- 
metry, simultaneously with taking a further step in estab- 
lishing the connection relations in non-Gaussian random 
matrix theory. Previously, the issue of parametric cor- 
relations in non-Gaussian random matrix ensembles (on 
the scale of several mean level spacings) has been consid- 
ered only by Hackenbroich and Weidenmiiller |l^ who 
used a supersymmetry formalism to provide a heuristic 
proof that in all three symmetry classes {(3 = 1,2 and 
4) the parametric spectral correlations of all orders in 
invariant matrix model given by Eq. (|l|) below are insen- 
sitive to the particular form of the confinement potential 
w[Ho], and thus follow the predictions found for Gaussian 
ensembles. This statement is valid for strong level con- 
finement and refers to the bulk of the spectrum. Until 
now, however, no constructive formalism was proposed 
to quantitatively describe parametric level statistics in 
a random matrix ensemble with a general non-Gaussian 
measure. Our aim here is to formulate an exact solution 
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to the problem of parametric level statistics for f3 = 2. 
Our study differs from Refs. |16 1^ in three main re- 
spects: (i) It applies to both strong and soft level con- 
finement; (ii) it is valid in arbitrary spectrum range; (iii) 
it gives a prescription of how to actually compute para- 
metric statistics for both arbitrary matrix size N and 
in the thermodynamic limit. As an illustration, we ap- 
ply the formalism developed to obtain closed analytical 
expressions for parametric eigenlevel correlations in the 
model Eq. for strong level confinement in the origin, 
the bulk, and the soft-edge scaling limits [Q. 

Now we turn to the formalization of the problem and 
to the statement of the main results. Consider an unper- 
turbed matrix ensemble defined by the joint distribution 
function of the entries oi N x N Hermitian matrix Hq, 



Pi[Ho]=Z^iexp{-Tri;[Ho]}. 



(1) 



Here v [Hq] is an arbitrary (not necessarily strong) con- 
finement potential ensuring existence of the partition 
function Z^. Eigenlevel correlations of arbitrary or- 
der in this ensemble are known j^] to be expressed 
through the scalar two-point kernel Kjsi {E, E') = 
^n=o (^) i^')- Here the fictitious wave functions 
i/)„ (£;) — A/'n exp {— t; (E) /2}p„ (E) involve a set of 
polynomials Pn (E) orthogonal on the entire real axis R 
with respect to the measure dii{E) — exp {—v (E)} dE, 



dlJL{E)pr, {E)pm {E)^Nn5. 



n^nm • 



(2) 



In the presence of the perturbation, which we param- 
eterize by the strength 0, the Hamiltonian of the sys- 
tem modeled by the matrix Hq flows to — Hq + 
(/)Hi. We choose the perturbation matrix Hi to belong 
to the Gaussian Unitary Ensemble (GUE), (57^ [Hi] oc 
exp I — Tr H^ } , and assume (f> to be positive without any 
loss of generality. Under these assumptions, the con- 
nected part of the parametric level densities correlator 
defined as 

{E, E'-^) = (Tr 5{E-li^)Tv5 {E' - H^)) 

- (Tr<5(^-Ho))(Tr<5(^'-H^)), (3) 

where the angular brackets (...) stand for averaging over 
the ensembles of random matrices Hq and Hi, is given 
by 
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^[v(E' +i4>ai)-v(E' +,j,a2)] 



Kn{E, E' + i(j)ai) 



X [Kn {E' + (t>(y2,E)-S{E' + 4>(y2-E)]. (4) 

Equation (0) is exact for arbitrary finite N , and the sum 
rule for R}ff is fulfilled by construction. Note, that for 
finite N , the equation contains both the two-point ker- 
nel Kn {E, E') of an unperturbed matrix ensemble Eq. 



(|l|) and the confinement potential v{E). It is worth men- 
tioning that only for the GUE measure, v\H.q\ oc Hq, 
this finite- solution reduces to the one known from the 
Dyson Brownian motion model |p8| in the limit (/) <C 1. 

Thermodynamic limit iV ^ cxd of Eq. (^) is substan- 
tially different in the case of the soft and the strong level 
confinement, (i) If the confinement potential is soft, the 
two-point kernel Kn has a well defined large-iV limit 
given by a corresponding Poisson kernel Kp{E, E';r) = 
Yln=o^''^n{E)ipn{E') taken at r = 1 As a conse- 

quence, Eq. (^) holds in the large-A^ limit as it stands, 
with Km replaced by Kp{r =1). After implementation 
of A^-independent unfolding, E ^ s, (j) ^ X, it provides 
an exact answer. Hence, generically there are no reasons 
to expect a universality of parametric level statistics in 
ensembles with soft level confinement, (ii) For strong 
level confinement, Eq. (^) can take a universal form in 
the thermodynamic limit after the implementation of the 
unfolding procedure ds/dE = Vj^, X — (pVjf (with Vpj 
being the local level density pG]), 



i?(^)(s,s';X) = -- 



da I 



da2 e 



xFx (s', ai,a2) K (s, s' + iXai) 

x[K{s' +Xa2,s)-6{s' + Xa2~s)], (5) 



where 



Fx is,ai,a2)= lim 



exp{^[v{EN,s + i'pN.xO'i)]} 



JV^oo expWviEN^s + (l)N,xcr2)]} 



(6) 



Equations (^ and (^ establish a link between the para- 
metric level statistics in unitary invariant non-Gaussian 
ensemble of random matrices and the energy level statis- 
tics in the same, unperturbed random matrix ensemble, 
through the unfolded two-point kernel 



K{s,s')= lim Km{Em,s,En.s') 

N—^oo 



xv^^''^{En.,s)vj^^''^{En.s') 



(7) 



Connection relations Eqs. (|j) - (0) provide a construc- 
tive rule for computing parametric level statistics in non- 
Gaussian matrix models and constitute the main results 
of the paper. 

In particular, when both Fx and K (s, s') are universal 
functions for a class of matrix models, the universality of 
eigenlevel correlations will automatically imply the uni- 
versality of parametric level correlations. To demonstrate 
this point, consider the matrix model Eq. (|l]) with strong 
confinement potential exhibiting a logarithmic singular- 
ity, v{E) v{E)-2a\og\E\ {a > -1/2). For /3 = 2 and 
in the absence of the perturbation, this model has been 
proven to possess three different types of locally universal 
eigenlevel correlations in the bulk (BSL) the origin 
(OSL) ||,| and the soft-edge (SSL) scaling limits j|. 
In all these scaling limits, the two-point kernel obeys the 
universal laws which we convert to 
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«'l(s)*2(s')-*2(s)^'l(s') 



(8) 



Here s and s' are appropriately unfolded spectral vari- 
ables, 

{TT^-^ sin (tts) , BSL, 
-(5/27r)i/2j„_i/2(7rs), OSL, (9) 
Ai (s) , SSL, 

1-2 (s) = *;(s) for BSL and SSL, while for OSL (s) = 
^'']^(s) — as~^^'i(s). These three universal kernels are 
referred to as the sine, the Bessel, and the Airy ker- 
nels, respectively. Universal parametric level statistics 
in this ensemble is a consequence of the connection re- 
lations Eqs. (^) and (||), which enable us to write down 
closed analytical expressions for the parametric level den- 
sity correlator i?^^^ just by noting that Fx is a univer- 
sal function: Fx{s, (Ji, (72) = (s + Xa2)°'{s + iXai) " in 
the origin scaling limit, and Fx = 1 otherwise. 

Let us now outline the derivation of the an- 
nounced relations. The two-point parametric correlator 
Rn {E, E'; (j)) of level densities can be represented as an 
average over the joint distribution function 



P2[Ho,H^] aexp{-Tr(w [Hq] 



H 2 



•P) )} ' 
(10) 



of two Hamiltonians, unperturbed Ho and perturbed H^, 
Rn [E, E'-^) = (Tr 5 [E - Hq) Tr 6 [E' - H^))^^ . (11) 

In Eq. (^1]) the averaging over the distribution 7^2 can 
explicitly be performed by using the Eynard-Mehta the- 
orem The theorem states that Rx (E, E'; (p) equals 
the determinant of the 2x2 matrix kernel. 



Rx {E,E';(j)) = Det 



a, 0=1,2 ' 



(12) 



where E^ = E + {E' - E) [I + (-1)"] /2. An explicit 
representation of the matrix kernel reads, 

C (^'^') = 'A"' E ^n'Pc.n (I) Q,,n (|) 



(13) 



Here Pa,n and Qa,n are four complete sets of orthogonal 
functions; two of them. Pi „ = P„ and Q2,n = Qn, are 
basic polynomials orthogonal with respect to the nonlocal 
weight W{i,ri) = exp {-u((/)^) - (C-??)^}, 



di / dr^Pn [0 W (C, v) (ry) = K5n 



(14) 



with hn being a normalization constant. (Hereafter, in- 
tegration over single variable runs over entire real axis). 



Two remaining orthogonal functions are connected to the 
basic polynomials as follows: 



P2,n (0 



dryP„(77)W^(ry,C), 



(15a) 
(15b) 



Remarkably, all of the four sets Pq.k and Qa.n can solely 
be expressed through the orthogonal polynomials p„ {E) , 
Eq. (^, which determine the eigenlevel statistics in 
an unperturbed random matrix ensemble Eq. (|lj). To 
prove this, we use the representation |Q of polynomials 
Pn and Qn orthogonal with respect to a nonlocal mea- 
sure as a determinant averaged over the distribution 
V2 [(/"Ml, 0M2], where Mi and M2 are n x n Hermitian 
matrices: 



PniO = Det(e-Mi) 



Qn (0 = Det (e - M2) 



7'2[0Mi,0M2] 



(16a) 
(16b) 



Integration over matrix M2 in Eq. (16a) is straightfor- 
ward; it yields the expression 

P„(0 = /Dct(e-Mi)\ ^cj^'/^pnm, (17) 

which we have recognized as a matrix-integral rep- 
resentation ]25[ | of the ordinary orthogonal polynomi- 
als appearing in E q. (p| ). To integrate out the ma- 
trix M2 in Eq. (16b) we rewrite the determinant 
there in the form of the integral Det„xn {£, — M2) = 
J D [x] exp {— {£, — M2) x} over n~component Grass- 
mann vector x — (xi; ■ • ■ j Xn)- Further integration over 
M2 leads to 



x(exp{-xt (C-Mi)x})p^[^M, 



(18) 



Decoupling the 'interaction' term exp{ — (x^x)^/4} in Eq. 
( |l8| ) followed by integration over the Grassmann vector 
X results in 



(0 



-1/2 



daie-^^P^i^ + iai) 



(19) 



Combining Eqs. ([l5|), ([l7| ) and ([l9|), we obtain explicit 
expressions for Pa,n and Qa,n entering Eq. ( p^ ) for the 



matrix kernel k 



(N) 



Pi.. {0 - 'i^'^'pn m , 

P2,„(e) =^'^^ J dvW{^,r^)pnm, 
-(^^)'/'e-'^(*«V„(^0, 



(20a) 
(20b) 
(20c) 



Q2,„(0 = / daie— Jp„(0(^ + iai)). (20d) 
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One can verify that the orthogonality relation Eq. jlj) 
is fulfilled, and the normalization constants in Eqs. ^) 
and (0 are related as /i„ = tt ^'^A/'„. Substituting Eqs. 
(^) into Eqs. (|l3| ) and (|lj), and subtracting discon- 
nected part given by the product k[f ^k2^^ is the final 
stage of the calculations which leads to our solution Eq. 
(^. Equations (||) and (^) follow upon appropriate un- 
folding. Note that these results can be extended to the 
n-point parametric correlators along the lines of Ref . . 

In conclusion, we presented an exact solution, Eq. (|^), 
to the problem of parametric level statistics in unitary 
invariant, non-Gaussian random matrix ensembles char- 
acterized by either soft or strong level confinement. The 
solution is exact in precisely the same sense as are the for- 
mulas by Gaudin and Mehta for conventional spectral 
statistics. Generically, the parametric 'density-density' 
correlator was shown to depend on both the two-point 
scalar kernel and the level confinement, through a double 
integral transformation which, in turn, provides a con- 
structive tool for description of parametric level correla- 
tions in non-Gaussian random matrix theory. In random 
matrix ensembles with strong level confinement, the so- 
lution presented takes a particular simple form in the 
thermodynamic limit leading to emergence of a new con- 
nection relation, Eq. (^, between the parametric spec- 
tral statistics and the scalar two-point kernel of an un- 
perturbed ensemble. In the case of soft level confine- 
ment, the formalism developed is potentially applicable 
to a study of parametric level statistics in unitary invari- 
ant random matrix model with log-squared confinement 
potential v{E) oc log^ \E\ at \E\ > 1. The latter 
problem is of conceptual interest in view of the recent 
work 101 where existence of a new class of random ma- 
trix ensembles with finite eigenlevel compressibility and 
multifractal eigenvectors was suggested. There, identi- 
cal two-level correlation functions in proper regimes in 
three different ensembles - (a) generalized matrix model 
of Ref. ||2^ , (b) power-law banded matrices of Ref. , 
and (c) rotationally invariant random matrices with log- 
squared confinement potential already mentioned above 
- were one of the key points to identifying the new uni- 
versality class relevant to the one expected at the edge of 
the metal-insulator transition. In this context, it would 
be desirable to learn whether these three matrix models 
also enjoy the same parametric level statistics. At least 
for two of the ensembles [(b) and (c)] such a verification 
seems now to be possible. 

I thank V. E. Kravtsov, K. A. Muttalib and I. Yurke- 
vich for discussions at different stages of this study. 
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